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Abstract 

We examine the order of the phase transition in the Potts model by using the 
graph representation for the partition function, which allows treating a non-integer 
number of Potts states. The order of transition is determined by the analysis of the 
shape of the graph-weight probability distribution. The approach is illustrated on 
special cases of the one-dimensional Potts model with long-range interactions and 
on its mean-field limit. 
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1 Introduction 



One of the less trivial questions in the study of phase-transitions is the de- 
termination of the order of transition. A large amount of work, especially on 
the Potts [1] model, ranging from the exact solutions [2,3] to different kinds 
of Monte Carlo (MC) simulations [4-12] (to cite only some of them) , has been 
done in order to establish reliable and at the same time applicable criteria for 
distinction between the first- and second-order phase transition. 

The coexistence of two different phases at the first-order phase transition 
point and its absence at the second-order phase transition point is the main 
physical fact on which all of the above methods (including the present one) 
are based. Performing the MC simulations on finite models in order to obtain 
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the probability distribution of some quantity (such as the energy or the order 
parameter) which has different values in different phases, is the usual way to 
investigate these phase transitions numerically [f 3]. In the present paper we 
show that the quantity which we define as the graph weight also has different 
values in the coexisting phases at the first-order transition point. In Section 3 
it is shown how the graph weight is related to the energy and the free energy 
of the system. The coexistence of phases with different energy at the first- 
order transition point then implies the coexistence of graphs with different 
weights. It can be identified by the two peaks in the graph-weight probability 
distribution. 

In order to test the graph approach, we analyse the Potts model in two cases: 
the mean-field (MF) case defined by the interactions of equal strength between 
all particles of the model and the one-dimensional (Id) model with long-range 
(LR) interactions decaying with distance r as r~^ 1+cr \ a > 0. In both cases two 
regimes (of the first- and second-order transition) are present. The regimes are 
separated by a point q c in the MF case and by a line (q c , a c ) in the LR case. 

The present approach has also the advantage of dealing directly with the 
non-integer values of the Potts states q. (Note that in reference [8] the non- 
integer q values were not obtained by a direct calculation, but by a histogram- 
like extrapolation from the results for integer-g.) This may be of interest in 
studying the threshold q c , separating the first- from the second-order transition 
regime, which is not always an integer. For example, in three-dimensional 
short-range (SR) Potts model, q c is not an integer [8,14], and in (lei) Potts 
model with LR interaction, the onset of the first order transition at q c seems 
to depend continuously on the interaction-range parameter a [15,16]. 

The plan of the paper is as follows. In the next section we define the model and 
introduce the corresponding graph expansion. The basic steps of MC algorithm 
are also explained in that section. In the third section we discuss separately 
the results for the MF and the LR interaction case. In the last section we 
summarize and discuss the advantages and open problems connected with the 
graph approach. 

2 Model and method 

The graph representation [17] was already used in MC investigation of critical 
behaviour of Potts models [5,6,10]. In the present approach we follow the 
line of reasoning given in reference [5], but with the basic difference that we 
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study directly the probability distribution of graph weights and not the cluster 
probability distribution. 

We begin by rewriting the model in the graph language. The reduced Hamil- 
tonian of the model, with periodic boundary conditions, has the form 

— H N ~ l N ~ i 

T 7f = E E Kj 5{ Si ,S i+j ). (1) 
K B -L i=1 j=1 

where Sj = 0, 1, . . . , q — 1 is the Potts particle placed on the i-th site of the 
chain, while 5 denotes the Kronecker symbol. K/s denote interactions between 
two particles at distance j. Due to boundary conditions in the LR case, each 
Kj has two contributions involving interactions at distances j and (N — j). In 
the MF case, the interactions are unique and the boundary conditions have no 
meaning. We use the system of units where Kj is proportional to the inverse 
temperature. 

The usual substitution, exp[Kj S(si, s i+ j)] = 1+Vj S(si, s i+ j), with Vj = exp(Kj) — 
1, leads to the partition function of the form 

N 9-1 N q-1 N-lN-i 

zn=u e exp(-ff/fcflT) = n e n n (2) 

1=1 s i= 1=1 Sl=0 1=1 j=l 

It is straightforward to establish one-to-one correspondence between each 
member of the r.h.s. of the above equation and a graph on the chain consisting 
of iV particles. Each square bracket can contribute to the above product in two 
ways: it gives 1 when there is no connection between the i-th and the i + j-th 
particle (inactive link), and Vj S(si, Sj +J ) when there is a direct connection (ac- 
tive link) between these two particles. The analytical expression attached to 
each graph Q will be called the graph weight, Wn(G)- The explicit expression 
for Wn(Q) follows from the r.h.s. of eq. (2). We can describe every graph con- 
sisting of c(G) clusters. By cluster we mean a set of particles interconnected by 
any type of active links and disconnected from other particles. Single particles 
are considered as one-particle clusters. The products of 5-functions from the 
expression for each graph will delete all except c g-summations on the r.h.s. 
of (2), so that the analytical expression for W N (Q) is 

W N {Q) = v b ^ G) v b ^ G) . . . v b N N :f ] q< g \ (3) 
The symbols bj(Q) denote the total number of active links of type j. 
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In this way, the summation over all 2 N ( N possible graphs between N 
particles corresponds to the partition function of the model 

Z N =J2 W N (G) = J2^n,w W N . (4) 

all g w 



The Mn,w denotes the number of different graphs with the same weight Wn- 
The above equation is thus the graph analogy of the more familiar expression 

Z N = J2^N, E e- E ^ T , (5) 

E 



where the summation runs over all different energies of the system, while 
Mn,e is their degeneracy. Indeed, behind the analogy, there is a connection 
between the average number of clusters and the active links and average of 
energy and free energy of the system. The derivative of the partition function, 
given by eqs. (3), (4) and (5), over temperature [10] leads to 

(E N )=T^~ 1 ^ 1 (b S )^, (6) 



while the derivative over q leads to 

/ dlnAr NtE (q) \ 1 / dE N \ 1 
\ dq T \ dq q 



(c). (7) 



The coexistence of two phases characterized by different values of energy cor- 
responds, by the above relations, to the coexistence of graphs with different 
weights. 

According to (4), one introduces the graph- weight probability distribution Pjv 



Numerically, Pn is obtained by a simple MC simulation of Metropolis type. 
The basic steps are: 

(a) Pick at random one link in the graph Q and change its status from active 
to inactive or vice versa. The resulting graph is called Q' . 

(b) Compare the random number < r < 1 with the ratio Wn(G')/Wn(G)- 
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If the ratio is smaller than r keep Q, otherwise save Q' as Q. 
(c) Return to (a). 

At each step, the number of clusters has to be counted. Since we deal with a 
model with interactions of infinite range, the determination of clusters could 
not be performed by determination of their boundaries, but one has to check 
for possible connections with all the particles of the considered cluster. To do 
so, during the simulations, we have to keep record of the cluster structure of 
the system, i.e. which particle belongs to which cluster. We start with a simple 
and known cluster configuration. In each of the following MC steps, when one 
link is changed, we examine whether this has produced a change in cluster 
structure: if the link is deactivated - whether the corresponding cluster is split 
in two or not; if the link is activated - whether two clusters get connected or 
not. 

Counting graphs by their weights through a large number of steps gives the 
unnormalized graph-weight probability distribution. Unfortunately, the short- 
coming of the graph approach, is that it takes a much longer time to obtain 
the same precision of results compared to the simulation techniques that can 
be applied only for integer q (like algorithm by Luijten and Blote [11]). On the 
other hand, the comparison of the present approach with the simple Metropo- 
lis single spin flip algorithm done for the model of N particles with integer q 
shows that 10 4 flips per link of everyone of N(N — l)/2 links in graph sim- 
ulations give approximately the same precision as 10 6 one-particle flips per 
particle in the spin simulations. For the chain of 400 sites considered here, 
the graph approach requires about 2 — 5 times larger CPU time than the 
Metropolis algorithm on spins. 



3 Results 



The simulations were applied to the two mentioned cases of the Potts model 
with N = 100, 150, 200, 250, 300, 350 and 400 particles. In final extrapolations 
only the data for N > 200 were used. The precision of simulations is deter- 
mined by performing 10 4 flips per link. 
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3.1 Mean-field case 



The MF case of Hamiltonian (1) is given by taking the equal strength of 
interactions among all the particles, i.e. Kj = K / N, where K denotes the 
inverse temperature. The exact work of Kihara et al. [2] puts in evidence that 
the MF Potts model has a second-order phase transition for q < q c — 2 and 
a first-order phase transition for q > 2. It thus qualifies as a good example 
for testing the present approach. In our simulation we start with the case 
q — 3. In continuation, we consider lower values of q, closer to the threshold 
q c in order to examine how efficient our approach can be in detecting weak 
first-order transitions. 

The shape of the finite- iV graph-weight probability distribution P N depends 
on temperature in a similar way as the energy-probability distribution does. 
Far from the transition temperature, Pn has a gaussian form. At the transi- 
tion, the shape of Pn depends on the order of the transition. For second-order 
transitions, the shape has a non-gaussian form, but it still has only one maxi- 
mum. For first-order transitions, Pn has two maxima which transform into two 
^-functions in the thermodynamic limit. We define the temperatures T/v where 
the two maxima in Pn are of equal height. The corresponding positions on the 
graph-weight axis are W N (for the ordered phase) and W N ° (for the disordered 
phase). When the two-peak structure of Pn becomes more pronounced with 
increasing N, one may conclude that the transition is of the first order in the 
thermodynamic limit with two coexisting phases: the one described by the 
graph weight W N — > W°, stable below the transition temperature T N — > T t , 
and the other one described by W N ° — ► W do stable above T t . 

In Figure 1 are presented the results of simulations for q — 3. It shows de- 
pendence of Pn on wn = In Wn / N In N at temperatures Tn- The peaks 
emerge for 200 < N < 400 and become more pronounced with increasing N. 
By the above arguments, the iV-dependence of P N confirms the existence of 
the first-order phase transition, in agreement with the exact solution [2]. 

The behaviour for smaller N considered (N < 200) points out that only one 
maximum in Pn may have different origins: either the transition is of the sec- 
ond order in the thermodynamic limit and we observe its finite- iV behaviour, 
or the transition is of the first order in the thermodynamic limit, but the cor- 
relation length, although finite, is comparable with the system size used and 
we are not able to see the coexistence of the two phases. Arbitrarily close to 
the threshold which separates the first- from the second-order transitions, the 
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correlation length is expected to become arbitrarily large and finally whatever 
system size is used, it will always be too small to show the first-order char- 
acter of the transition. Such a situation is expected in the MF case when q 
approaches 2. The simulations on q = 2.8 model performed at the correspond- 
ing temperatures defined as Tjv give the distributions similar to those shown 
in Fig. I. On the contrary, simulations performed for q — 2.1 and 2.5 do not 
show the two-peak structure in the graph-weight probability distribution for 
the considered system sizes. 

To some extent, the above presented simulations can be compared to the exact 
solution. Since A/at,e(5) and E^{q) can be exactly calculated for the transition 
temperature T^ F = 2[(q — l)/(q — 2)] In (q — 1), [2], it is easy to obtain from 
eqs. (6) and (7) the analytic expressions for the positions ( w^ do ) of the two 
peaks at temperature Tmf in the large N limit: 



N 



, M9 -^r M '- 1)+ i^v +c V' 

"~ 1 ln(9-D+ a « 



do 



q(q-2) 



+ o(— 

In N K N' 



(9) 



where 



( q - iy + 1 



fq 5 2 (q - 2) 

do q- 1 



ln(q — 1) In 



q (q ~ 2) 



ln(g — 1) In 



2 q —± ln(q - 1) 

q z 

q-l 



+ 



\nq 



Hq-1) 



q-2 



\n(q- 1) 



+ In q 



5(5-2) 
(q — 1) ln(g — 1 

5(5-2) 



In Table 1 we compare the data for q = 3 obtained by MC simulations at 
T/v =4 oo = 0.3566, with those calculated from (9) for finite iV = 400 and for 
N -> oo. 
Table 1 



Simulated (MC) and exact values of the positions ( w ° N ) and ( w do 



N 



of two maxima 



in graph-weight probability distribution of q = 3 MF Potts model. 



T 


N 


Method 




«> 


T N 


400 


MC 


-0.5271 


-0.3123 


Tmf 


400 


exact 


-0.5045 


-0.2847 


Tmf 


oo 


exact 


-0.6931 


-0.4621 



The values of ( w^ do ) cited in Table 1 show that the simulated results fit well 
the exact values for finite N = 400. Relatively high difference (24% for (w^ ) 
and 32% for {wfj>)) compared to the exact iV — > oo values comes from the 
slow 1/ln N convergence. 
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Table 2 

Values of MF temperatures Tjy F and LR temperatures T^ R at which P/v shows two 
peaks of approximately same heights, followed by the extrapolated values. 



N 


/tiMF 
1 N,q=2.S 


1 N,q=3 


1 N,q=3,<r=0.1 


200 


0.3694 


0.3533 


3.25248 


250 


0.37042 


0.35491 


3.3402 


300 


0.3699 


0.35490 


3.4235 


350 


0.37298 


0.35615 


3.4880 


400 


0.3730 


0.3566 


3.55304 


oo 


0.376 


0.360 


7.73 



In Table 2 we present the corresponding data for T/v (obtained in simulations) 
in function of size. Since for large N the dependence of TVs is approximately 
linear in iV -1 , we use the extrapolation form T N = T M F,t+b/N. The coefficients 
T MF ,t and b were obtained in least-squares approximation (LSA). By taking 
into account only the data for sizes 200 < N < 400, one obtains T M F,t(<? = 
2.8) = 0.376 and T MF j(q = 3.0) = 0.360. In both cases, the difference between 
the extrapolated values and the exact ones T MF (q = 2.8) ~ 0.378 and T MF (q = 
3.0) ~ 0.361, is less than 1%. 

3.2 One- dimensional long-range case 

In the case of one-dimensional Potts model with interactions decaying as 
l/r 1+a , the interaction strength Kj is given by K [1 / j 1+a + 1 / (TV - j) 1+a ], 
where K denotes the inverse temperature. This model has a phase transition 
at finite temperature [18,19] for all q and < a < 1. Our recent MC simula- 
tions of the energy probability distribution on q = 3 and q = 5 models [15,16] 
have shown that the order of the transition depends on q and a. For the same 
fixed value q — 3, it was shown that the transition changes from the first- to 
the second-order one with increasing a. 

Again, according to the relations (6) and (7), one can interpret the meaning 
of the graphs through their connection with the energy and the free energy of 
the system. 

The purpose of simulations presented in this subsection is twofold: first, we 
wish to test the above described graph approach comparing the results with 
those obtained earlier by MC simulation of energy probability distribution 
(strong first-order transition for model with q = 3 and a = 0.1); second, by 
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using the same graph approach we wish to analyze an example of a non-integer 
q model with a second-order phase transition [20,21] (we choose q = 0.5 and 
a = 0.1 and 0.8). 

The result of simulations for the case with q — 3, a — 0.1 is shown in Fig. 2. 
It presents P N versus w N for system sizes ranging from 200 to 400, at tem- 
peratures T N where two peaks of the approximately same height appear. The 
figure shows that the depth of the minima increases with N . According to 
the discussion in the preceding subsection, a behaviour like that is charac- 
teristic for a first-order phase transition. This conclusion is confirmed by the 
results[16] of the MC simulations for the energy probability distribution. 

The temperatures T/v, where 200 < N < 400 are presented in Table 2. Com- 
pared to the MF results and also to those for higher values of a, the results 
for a — 0.1 converge very slowly, so that, for sizes considered here, we do not 
expect so good accuracy in extrapolation to N — > oo. The extrapolation which 
takes one correction term, of the form T N = T LR:t + b/N x gives T LRyt = 7.73 
with the convergence exponent x = 0.1, while taking the form with two cor- 
rection terms, e.g. T N = T LRjt + b/N x + c/N 2 ' x gives by LSA fit x = 0.12 
and Tl R j = 7.45, value which differs by 4%. The first result for Tlrj shows 
the discrepancy of 8% compared to the improved finite-range scaling (FRS) 
[19] result (T FRS = 7.14) cited in [16], and even larger discrepancies when 
compared to the renormalisation-group (RG) result (T RG = 6.72) [22] and the 
earlier MC result (T MC = 6.25) [16]. (Notice however that, due to the same 
slow convergence, the difference between these previous results is also quite 
significant when a = 0.1.) Apart from the reduced precision due to the small 
convergence exponent, the additionnal source of deviation can come from the 
crossover effect and be the consequence of the limitation to relatively small 
sizes for which the system does not exhibit full qualities of the first-order 
transition. 

In order to investigate the case with non-integer q, we also considered the 
case q = 0.5 with a = 0.1 and 0.8, where by [20] and [21] the transitions are 
of classical and non-classical second-order type, respectively. The simulations 
were performed for a wide range of temperatures around critical temperatures 
obtained earlier by FRS [19]. As one can expect on the grounds of the discus- 
sion at the beginning of the preceding subsection, no double-peak structure in 
graph-weight probability distribution has been found. 
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4 Conclusion 



The graph-weight probability distribution is introduced and applied to the 
analysis of the order of transition on two special cases of the Potts model: 
its mean-field case and the case in Id with power-law decaying interactions. 
The mean-field case was analyzed for q = 3 and several non-integer values 
approaching q c = 2. The long-range case was examined for q = 3 and q = 0.5 
with o = 0.1 and o = 0.1,0.8, respectively. The simulations were limited by 
time to the systems of sizes N < 400. 

It is shown that the graph weight is an appropriate quantity to distinguish the 
coexisting phases at the first-order transition point. The physical interpreta- 
tion of the graph weight becomes transparent through the equations (6) and 
(7) which relate the average number of active links and clusters to the average 
of energy and the free energy of the system. 

By analyzing the graph-weight probability distributions of the MF and Id LR 
Potts models, we have observed the behaviour of distributions characteristic 
for first-order transitions in the case of q = 3 and 2.8 MF models and q = 
3, a = 0.1 LR model, while the transition of the second order was obtained in 
the case of q = 0.5, a — 0.1, 0.8 LR model. 

Transition temperatures have also been calculated. The estimated values of 
those temperatures in the thermodynamic limit agree for the MF case values 
within 1%, while for the case with power-law decaying interactions, where only 
the approximate results are available, it agrees with the discrepancy of 5% to 
24% depending on the method. 

All of the above facts qualify the graph-weight probability distribution as 
an alternative quantity in investigations of the order of transitions in Potts 
models, capable to deal directly with non-integer values of q. This makes the 
graph approach interesting in the framework of efforts of determination of the 
border between the first- and second-order regions by continuously varying q 
in the (q, d) plane of SR models or (q, a) plane in Id LR models. 
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Figure captions 



Fig. 1: The simulation of the graph- weight probability distribution P/v of 
the q = 3 MF Potts model versus = In Wn / N In N performed at 
temperatures T/v (see in text). 

Fig. 2: The simulation of graph-weight probability distribution of the 
q — 3, a — 0.1 Id LR Potts model versus = In W 7 ^ / N In iV performed 
at temperatures T N (see in text). 
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Fig. 1 (q = 3, M F) 
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Fig. 2 (q = 3, a = 0.1) 
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